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1 Introduction 



Lattice Monte Carlo simulations provide a powerful tool in studying non-perturbative effects in 
QCD. Promising applications are, for instance, hadron spectroscopy, high temperature thermody- 
namics or the structure of the QCD vacuum. However, using a space-time lattice as an ultraviolet 
regulator introduces discretization errors. A careful analysis is required to extrapolate the Monte 
Carlo results to zero lattice spacing recovering continuum physics. The simplest choice of the 
lattice action corresponding to the given continuum theory produces, in general, quite large lat- 
tice artifacts. It is advisable to introduce more complicated lattice actions which can significantly 
reduce such discretization errors. 

Wilson's renormalization group (RG) approach jl], 1 1 offers a promising method for eliminating 
lattice artifacts in asymptotically free field theories Using actions on the renormalized 

trajectory (RT) of a given RG transformation (RGT) one eliminates lattice artifacts in physical 
quantities. These were called perfect actions The RT starts from a fixed point (FP) on the 
critical surface. The FP action is defined by retaining the functional form of the FP and using it 
also for finite couplings. In other words, the FP is a point in the space of couplings (including the 
prefactor (3 oc l/g 2 ) while the FP action is a straight line going in the direction of g 2 ). The FP 
action is a good approximation to the RT. It has the important property that it defines a perfect 
classical lattice theory: The solutions to the lattice equations of motion have exactly the same 
value of the lattice action as the continuum action on the corresponding continuum solutions. For 
example, the FP as a classical action has scale invariant instanton solutions when the corresponding 
continuum theory has. Furthermore, it can be argued that the RT coincides with the FP action 
even in 1-loop perturbation theory — the FP action is '1-loop perfect' |^,|fj,[|. In the last years 
the FP action has been constructed for several two- and four dimensional asymptotically free 
theories §-Q. 

The type of blocking used in the RGT is to a large extent arbitrary. But the properties of 
the FP, most notably the interaction range, are very sensitive to the choice of the parameters 
in the blocking transformation. It has been demonstrated that one can choose the RGT in such 
a way that the FP action is essentially concentrated on the unit hyper-cube. This is important 
since one has to simplify the FP action in order to use it in numerical simulations. In measuring 
physical quantities, like masses, the string tension or critical temperature approximate FP actions 
seems to give much reduced cut-off effects. But in expectation values of local operators — like 
the static gg-potential — new lattice artifacts appear. To reduce such cut-off effects one has to 
use for a given FP action (or even for an action on the RT) improved local operators. In the RG 
approach the origin of these artifacts is quite obvious. Although the partition function and all 
long range properties (masses, etc.) are unchanged by a RG step, the fields on the coarse lattice 
represent some averages over the original fine fields. If the blocking procedure violates rotational 
invariance, correlation functions of the coarse fields will show a violation of this symmetry as well. 
(As an analogy, in continuum electrodynamics the potential between two square-shaped charge 
distributions is not given by the 1/r potential alone, it has contribution due to higher multipole 
moments as well. Still the physical content is described exactly, only the testing objects are chosen 
inconveniently.) One way out, as mentioned above, is to construct better fields (FP fields) - 
another is to use a better, i.e. more rotationally symmetric RGT to make the testing objects, 
the bare fields, more spherical. The most important requirement remains, of course, the short 
interaction range of the corresponding FP action. 

The present work has been motivated by the fact that the pcrturbative potential of the bare 
Polyakov loops for the FP actions considered in ref. |], (jj has shown a considerable violation of 
rotational symmetry at a distance of 1-2 lattice units. Our aim is to extend the class of blocking 
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transformations to contain 'more rotationally invariant' ones. By optimizing the free parameters 
of the transformation for a short interaction range, we obtain a FP action which has, as expected, 
also smaller violation of rotational symmetry in the correlators of bare Polyakov loops. Apart from 
the complications connected with the more general form of the RGT this work repeats the steps 
described in detail in ref. |],[7| and we shall try to avoid unnecessary repetitions. 

While a parametrization of the FP action should represent the true FP action sufficiently well, 
it should also be simple enough to be used in numerical simulations. Parametrization errors are 
one source of 'imperfectness' in the present approach. Another is that the FP action is only an 
approximation — although a very good one — to the perfect action, i.e. an action on the RT. The 
latter can be improved by trying to follow the RT |j|,[l3],[l7]], the former by including more operators 
in the parametrization. How far one should go in this direction depends on many things — how 
large are the remaining lattice artifacts, what is the overhead for including further operators, what 
is the gain in computer time by working on coarser lattices, etc. A detailed study of the cut-off 
effects is needed to decide how a given lattice action performs. In this paper only the first few 
steps are taken in this direction for the new RGT. 

The Wilson action in the SU(3) gauge theory gives huge cut-off effects for thermodynamic 
quantities (e.g. about 50% for Nt = 4, (3 — ► oo), therefore it is especially important to use 
improved actions for such studies |lj,^9). In a related paper A. Papa has studied lattice artifacts 
in thermodynamics with the new FP action considered here p0|] . 



An alternative way of improving lattice actions is the Symanzik approach J2 1| — 124[| . Here one adds 
to the action order by order in perturbation theory appropriate irrelevant operators which cancel 
the leading cut-off effects. Although it is based on perturbation theory, the idea can be extended 
beyond the perturbative region as well |25| ]. The tadpole improvement also aims at extending the 
Symanzik approach beyond the perturbative regime, using a non-perturbative phenomenological 
determination of the coefficients of the correction terms p6|,|27|. 

We note here that the use of FP actions is especially promising in the study of topological 
effects |I^,|lJ,^6[^8| , because in this approach there are no dangerous 'dislocations' and an improved 
definition of the topological charge can be given. 

The paper is organized as follows: First we review the RGT and its FP with special consid- 
eration to the new block transformation. We study the quadratic expansion and determine the 
optimal block transformation parameters. This is followed by the construction of three different 
loop parametrizations of the FP action. In section || we study cut-off effects in the quadratic 
approximation of the static gg-potential comparing the Wilson action with the FP actions derived 
earlier and the new FP action. To set the physical scale we measure the critical couplings for 



various lattice sizes. This is done in section 3.2. We then study the performance of the new FP 



action by evaluating the static gg-potential using Polyakov loop correlators. Finally, in the last 
section we present some conclusions. 



2 The FP action 



With emphasis on the new block transformation we briefly review the RGT and its FP presented 
in ref. || . We evaluate the FP action in a quadratic approximation and fix the free parameters in 
the block transformation. Finally, we construct simple parametrizations of the FP action, one of 
which will be used in MC simulations. 
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We consider a SU(N) pure gauge theory]^] in four dimensional Euclidean space denned on a 
periodic lattice. The partition function is defined as follows 

Z= [Due-? 8 ™, (2.1) 



where DU is the invariant group measure and (3S(U) is some lattice regularization of the continuum 
action. Starting from such a lattice action we perform a real space RGT by defining a new action 
by 

e -p's'(v) = [ DUexp[-P{S(U)+T(U,V)}}. (2.2) 



Here T(U, V) is the blocking kernel of the transformation and is defined as 

T(U,V) = £ (X^b) - £ Re Tr [ ^(n B )Qt(n B ) ] ) , (2.3) 



N 

where Q^(tib) is a JV x iV matrix which represents some mean of products of link variables U^n) 
connecting the sites 2ns and 2ns +2/5 on the fine lattice. The term N^ifis) is a normalization 
which ensures the invariance of the partition function. The parameter k is a positive constant. It 
will be fixed later to optimize the quadratic lattice action. 

The block transformation is defined by a prescription for constructing the averaging matrices 
Quins)- Two such block transformations were defined in rcf. 01. Here we introduce one, which we 
call type III blocking. Instead of using just simple staples we also build 'diagonal staples' by going 
first in the planar or spatial diagonal directions orthogonal to p, then after a step in direction p 
returning along the corresponding diagonal. To be specific, first we create the matrices W^ m '(n, n') 
connecting the sites n and n', where n' is a site with coordinates n'—n^ — 0, \n' v — n v \ < 1 (for 
any v): 

W (Q) (n,n) = 1, (2.4a) 
W [1) {n,n + v) = U u {n), (2.4b) 

W {2) {n,n + v + p) = ^{U v {n)U p {n + u) + U p {n)U v {n + P )) , (2.4c) 

W^in, n + i> + /5 + A) = i (U u (n)U p (n + v)U\(n + v + p) + permutations) . (2.4d) 

Here v, p and A go over all (positive and negative) directions different from p and from each 
other. (Of course, U— V {n) = U\{n — v ). Values of W {m ^(n,n') not indicated are take n to b e zero. ) 



(n, n') represents the 'planar diagonal link', (n,n r ) the spatial one. In eqs. (2.4c), (2.4d) 
the sum is taken over all shortest paths leading to the endpoint n' of the corresponding diagonal. 
The fuzzy link operator is constructed then by a modified Swendsen smearing p9| : 

3 

W,(n) = ^ ^ Cm ^(n, n ')C/ M K)W (m) K + ^« + ^- (2-5) 

m— n' 

The coefficients c m are free parameters subject to the constraint: 

c + 6ci + 12c 2 + 8c 3 = 1. (2.6) 

Finally, the matrix Q^ns) is the product of two fuzzy link operators connecting the points 2ns 
and 2ns + 2p on the fine lattice: 

Qn(n B ) = W A1 (2n s )W p (2n s +/i). (2.7) 



1 Although the numerical studies have been performed for SU(3), the equations are mostly given for general N. 
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The condition eq. (2.6) ensures that for a trivial field configuration Q^(ub) is equal to the unit 
matrix. 



In the limit (3 — > oo eq. (2.2) reduces to a saddle point problem. At the FP this is an implicit 
equation for the FP action S FF : 



S FP (V) = mm L FP (U) + £ (^(n fl ) - £ Re Tr [V,(n B )Ql (n B )] ) ] 

I n B ,H ) 



(2.8) 



The normalization constant at j3 — oo becomes 



- Re Tr[^ Q t ( n s )]} . (2.9) 



As was already shown in ref. |3], p[ the FP action has scale invariant instanton solutions. This 



is a consequence of the FP equation fl2.8| ) alone and applies also for type III blocking. We will not 
proceed along these lines however, but will turn to an expansion in the vector potentials. 



2.1 The quadratic approximation 

Again this section is very similar to the one presented in ref. || . The only difference to type I and 
type II transformations is the actual transformation tensor involved. 



The FP equation (2.8) is valid for arbitrary configurations {V}. For smooth configurations it 
can be expanded in powers of the vector potentials B^ub) and A^(n): 

V„(n B ) = e iB ^ nB \ U^n) = e iA ^ n \ (2.10) 



The most general form of the quadratic action can be written in Fourier space as 



2NS(u) = ^J2p^ Tt M- k )Mk) 



V 



OA 



(2.11) 



where V is the volume of the fine lattice and p^vik) are the quadratic coefficients to be determined^. 
The transformation kernel becomes 

2NT(U, V) = V Tr [(T^-kg) - B^-kg)) {T^kg) - B^kg))} + O (cubic) , (2.12) 

VB , — 



where Vg = V/16 is the volume of the coarse lattice. The matrices T^kg) represent the linear 
contributions to the averages Q^{ns) and can be written as 



(2.13) 



1=0 



where I = (Iq, l±, I2, 1 3) is an integer vector and the summation goes over L = 0, 1. The tensor lu^ 
is fixed by the form of the blocking kernel: 



uJ^(k) = (1 + e ife ") \co5„ v + 6cirW (As) + 12c 2 t$ (k) + 8c 3 r<$ (fc) 



(2.14) 



In the following we suppress the tilde for Fourier transformed quantities. 
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Here the coefficients c m were introduced in eq. (2.5). The quantities T^\k) are the linear contri 



butions to the fuzzy link operator; they can be written as 



6 
1 

24 
1 

48 



KK(8 - £ + k^k* + k v k* v ) + $ M „(24 + ZKK) 



^(24 - + r - ?7 + 6fc^ - 2^* + 2(k„k*Y + Qk v k* v 
— 2£fc„fc* + 2(k„k*) 2 + 2fc^/c* k^k*) 



5^(48 - 24£ + 6£ 2 - £ a - 6£fc M fc; + 2f%*c* - 2£(^fc;) 2 
+ 3£r? - 2e - 6ry - 2i 7 fc /1 fc*; 



The lattice momentum fe„ is defined by 



kfj, — ' 1 



and we have introduced the following abbreviations 

iu n n 

As a consequence of gauge invariance rju?\k) and cj m „(/c) satisfy the relations 

t"u!/ {k)k v = A:^ , Ci^,„(fc)fc„ = 2fc^ . 



(2.15) 



(2.16) 
(2.17) 

(2.18) 



Combining the above expansions the FP equation has at the quadratic level the following form 
^T.P'^(k B ) TT[B^-k B )B u (k B )}^mm[^Y^p^(k) Tr [A„(-k)A„(k)] (2.19) 

+ ^ E Tr KM-M - (T M (fc B ) - B,{k B ))] }. 



To proceed further we have to introduce a temporary gauge fixing to be able to invert the 
tensors p^ and p'. One obtains then from eq. ( [2.19 ) a recursion relation for the propagator 



D„u{k) = p~l{k). 



(2.20) 



A way to find the FP solution for the propagator is to start from the Wilson propagator with some 
gauge fixing parameter a: 



\k\ 2 |fc| 4 ■ 



(2.21) 



and iterate the propagator to the FP. (At the end we have to switch off the gauge fixing in p^, 
which means taking the limit a — > oo.) Using gauge relations for (^^(fc) one can show that 
starting from the standard propagator, after an arbitrary number of iterations, the propagator 
Dfj,i,(k) assumes the form: 

D^{k) = G^(fc) + af{k)%k* v , (2.22) 
where G^(k) and f(k) are independent of a. Under the RGT they are iterated as follows: 



16 ^ 

1=0 



( ks , v\n( kB , n \( kB , i\ 



1 



(2.23) 
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r 




Poo(r-) 




r 




Pio(r) 











1 


-0.47938 














-0.72414 








1 


1 


-0.10005 











1 


-0.05969 





1 


1 


1 


-0.03110 








1 


1 


-0.00991 


1 








1 


0.00588 


2 


1 








-0.00222 


1 





1 


1 


-0.00457 





1 





1 


-0.00204 


1 


1 


1 


1 


-0.00316 


2 


1 





1 


-0.00081 


1 








2 


0.00208 





1 


1 


1 


-0.00059 


1 











0.00204 


2 


1 


1 


1 


-0.00043 


2 











0.00138 





2 








-0.00023 











2 


0.00132 











2 


0.00022 








1 


2 


-0.00084 








1 


2 


-0.00011 



Table 1 : Some of the elements of poo ( r ) and Pw (r) in configuration space for the RGT of type III. 
The values are obtained after 5 RG steps on a 16 4 lattice. 



and 



1=0 



For finite a the action density becomes 



af(k)u(k) 



where we have introduced 
In the limit a — > oo one obtains 

P^{k) = G-J,{k) 



u(k) = KG-hk)k p . 



G- l p {k)k p -KG-l{k) 
u{k) 



(2.24) 

(2.25) 
(2.26) 

(2.27) 



The type III RGT has four real parameters. These are ci, c%, C3 and k in eqs. (2.5) and 



(2.3). We tune them for a short ranged quadratic FP action. Since we optimize in a larger set 
of blocking transformations (note that at the quadratic level the blocking transformation includes 
those considered in type I blocking), the final action is expected to have a shorter interaction range. 
We found the following parameters to be optimal: 



ci = 0.07, c 2 = 0.016, c 3 = 0.008, k = 



(2.28) 



The iteration procedure converges quite rapidly with a next-to-leading eigenvalue of 0.25. While for 
type I and type II blocking the largest couplings of p p „{r) decrease as exp(— 2.7r), rsp. exp(— 3.1r), 
for type III blocking it is exp(— 3.4r). Some of the values of p^jf (r) have been listed in table [j]. 



Solving eq. (2.19) we obtain a linear relation between the field on the coarse lattice and the 
minimizing field on the fine lattice: 



where 



Mk) = Z MV (k)B v (2k), 
Z^(k) = [D(k)^(k)D'-\2k)] ^ 



(2.29) 
(2.30) 
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r 




Zoo{r) 






r 




Z w (r) 
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0.12507 





1 








-0.03993 








1 


1 


0.06061 


-1 


1 








-0.01336 


-1 











0.05378 


-1 











-0.01237 


-1 








1 


0.04776 





1 


1 





-0.01164 





1 


1 


1 


0.03794 





2 








-0.01053 


-1 





1 


1 


0.02692 





2 





1 


-0.00485 


-2 











0.02274 


-1 


1 





1 


-0.00484 


-1 


1 


1 


1 


0.01785 





1 





2 


-0.00451 











2 


0.01165 





1 


1 


1 


-0.00416 


-2 








1 


0.01058 


-1 


2 








-0.00409 


-1 








2 


-0.00963 











1 


-0.00406 



Table 2: Some of the elements of the tensor Z^{r) in configuration space. The lattice size is 16 4 
and 5 RG steps were taken. 



The connecting tensor Z^ v enters, in particular, in the construction of FP operators ||. For large 
a it is given by: 



Z, v {k) = [G(fcVt(fc) G '- 1 (2fc)]^ + — • 2kl [G'-\2k)] av 



oil 



(2.31) 



where u'(k) is defined in eq. (2.26) with G replaced by G' . The gauge relation for Z^ reads 

Z^ v (k)2k v = kfj,. (2.32) 
Some of the Z^(r) values are listed in table ||. 



2.2 Beyond the quadratic approximation 



For the use in numerical simulations we need a parametrization of the FP action S FP in terms of 
gauge invariant products of link variables (n) . For type I and type II RGT this has been studied 
in detail in ref. Q|. Here we shall simply outline the procedure and present the results. 



We parametrize the FP action by powers of traces of loops: 



(2.33) 



C m 



where Uc denotes the product of link variables U^n) along the closed path C. For a set of quadrat - 
ically independent loops the coefficients c\ (C) can be obtained from the quadratic approximation 
Pnv(f)- We apply a numerical procedure to determine those coefficients with m > 2. We generated 
about 500 configurations V using the Wilson action with j3 ranging from 5.1 to 50.0 and deter- 
mined the corresponding fine configurations U(V) by numerical minimization. The procedure is 
simplified by the observation that the typical value of the action density on the minimizing con- 
figuration U(V) in eq. (2.8) is by a factor 30-40 smaller than the action density on the coarse 
configuration V. This allows us two make a two step iteration: First by using a good quadratic 
approximation on the fine lattice one determines a parametrization of the l.h.s. which includes 



7 



higher powers and describes the 'measured' values of the action well (i.e. those obtained by the 
minimization procedure) . In the next step this precise intermediate parametrization is used on the 
r.h.s. of eq. fl2.8| ) to determine the value S FP (V) for any configuration V, including the very rough 
ones. The final step is to represent the set {V^, S FP (V 1 -^)} by some form suitable for numerical 



simulations. Considering a larger number of loops in eq. (J2.33) one gets a better parametrization, 



but the computer time grows rapidly with the length of the loops, hence one has to restrict the 
number of loops. 

Below we present three alternative parametrizations of the FP action of type III blocking. 
They are denoted as type Ilia, Illb, and IIIc. As has been mentioned earlier [|[[| there are many 
parametrizations which represent the given set of action values equally well. The optimal values 
are not unique — the matrix in the corresponding \ 2 fit has many nearly zero eigenvalues. Since 
the roughest configurations dominate the fit our results will represent the FP action on typical 
configurations generated with /?wiison ~ 5 — 6. 

We shall consider the three loops shown in fig. [j]. In the parametrizations Ilia and Illb we 





pl 



rt 



Pg 



Figure 1: Three different loops considered in the parametrization of the FP action. The plaquette 
operator (pl), the rectangle (rt) and the three-dimensional parallelogram (pg). 



have chosen the same set of operators as in ref. M for the type I FP action, the plaquette and 
the rectangle. For type IIIc action we have also included the rectangle. The maximal power m 
in eq. ( 2.33 ) for all loops was 4. For the type Illb and type IIIc parametrization we have used 
the freedom in the coefficients to demand that they satisfy the tree level Symanzik condition for 
spectral quantities |||]: 

ci(pl) + 20ci(rt)-4ci(pg) = 0. (2.34) 
The results are listed in table ||. As stated above all the three fits represent the FP values of our 





operator 


Cl 


C2 


C3 


c 4 


Ilia 


pl 


0.4822 


0.2288 


-0.1248 


0.02282 




Pg 


0.06473 


-0.02245 


0.002954 


0.003471 


Illb 


pl 


0.3333 


0.4845 


-0.2730 


0.04858 




Pg 


0.083333 


-0.01644 


-0.009525 


0.006324 


IIIc 


pl 


0.4792 


0.2226 


-0.1273 


0.02403 




rt 


-0.0091 


-0.04471 


0.02563 


-0.003698 




Pg 


0.0742 


0.02047 


-0.02398 


0.007730 



Table 3: Couplings of the 8 parameter FP action of type Ilia and type Illb, and the 12 parameter 
FP action of type IIIc. 

set of configurations equally well. The deviations are similar to those shown in ref. 0. We expect 
the three types of parametrizations Ilia, Illb and IIIc to give comparable results for simulations 
in the range /3wiison ~ 5 — 6. At small lattice spacings, however, when the linear approximation 
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dominates, the Symanzik condition imposed on type Illb and type IIIc actions will eliminate the 
0(a 2 ) effects, which makes them better suited for perturbative calculations. The couplings in 
table ^ also satisfy the normalization condition 

d (pi) + 8ci (rt) + 8 Cl (pg) = 1 . (2.35) 



Note that the whole quadratic part in the full FP action satisfies exactly the Symanzik con- 
dition. However, in a parametrization using a restricted set of loops this condition may not be 
satisfied any more — unless it is specially requested as an auxiliary condition. The quadratic 
part (the coefficients C\{C)) used in the parametrization IIIc are obtained from the best fit to the 
quadratic coefficients p^vir) — with those operators under the constraint ( 2.34 ). All the numer- 
ical simulations have been performed using the parametrization Ilia. The other two forms are 
given only to show the ambiguity of the truncation step and to provide a parametrization which 
maintains the Symanzik improved nature of the full FP action. 



3 The static gg-potential 



A good check for discretization errors is offered by the static gg-potential. Below the critical 
temperature T c the potential is expected to rise linearly for larger distances, with the string tension 
as its slope. For smaller distances the potential should exhibit a Coulomb potential. In lowest 
order perturbation theory it is proportional to the time average of the zero component of the free 
field propagator. In the following two sections we study the perturbative potential to lowest order 
and the full potential using Polyakov loop operators. 



3.1 The tree level potential 

At the quadratic level the static potential in a finite continuum boxQ of size L^. is proportional to 
the Coulomb potential: 

v cont (x) = ^ yy^-^' where k i = Y^k- (3.i) 

On a lattice of size Nf and an arbitrary lattice action we get 

V ^ = ^3 E e ^A)o(fco = 0, fc), (3.2) 

s fc#0 

where D^ v (k) is the propagator corresponding to the quadratic part of the lattice action. On fig. 
we plotted the differences between the lattice potentials and the continuum potential for Wilson 
action and type I, II, and III actions. To see how large this deviations are one can compare the 
force measured on the lattice to the continuum force. In table || the relative differences 

R{ri,r 2 ) = prr 77 prr (3.3 

VWit(n) - V con t(r 2 j 

are listed for the few first distances on the lattice. There are two types of deviations from the 



iNote that the force at \x\/L 1/3 in a finite box differs from the infinite volume result by as much as 25% so 
it is essential to use the correct finite box expression for comparison to the lattice results. 
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n 


T2 


Wilson 


type I 


type II 


type III 


(0,1,1) 


(0,0,1) 


0.325 


0.2797 


0.1479 


-0.0080 


(1,1,1) 


(0,1,1) 


0.125 


0.0051 


0.0113 


-0.0044 


(0,0,2) 


(1,1,1) 


-0.895 


-0.0544 


-0.0725 


0.0125 


(0,1,2) 


(0,0,2) 


0.662 


0.0034 


0.0493 


-0.0048 


(1,1,2) 


(0,1,2) 


0.283 


0.0092 


0.0118 


-0.0022 


(0,2,2) 


(1,1,2) 


-0.103 


-0.0016 


-0.0022 


0.0004 


(1,2,2) 


(0,2,2) 


0.183 


0.0017 


0.0023 


-0.0005 



Table 4: The relative error of the perturbative force on the lattice (cf. eq. (3J3)) for various lattice 
actions. 



continuum result for the FP action ||. The first type is due to the deviation of the shape of 
the blocking from a spherical averaging. This part contains, for example, a term oc P^{cos9) /r b 
(an octupole term) and can be corrected by using FP Polyakov loops^j. The second type is an 
exponentially falling correction cx exp(— r/r^) where ro is of the order of the interaction range and 
is due to the effect of the quantum fluctuations on the fine lattice. As seen from fig. 0, for type I 
and type II actions the discretization errors are much reduced, but at r — 1 a substantial deviation 
from the continuum can be observed. On the other hand, for type III action the perturbative 
potential lies for all distances very near to the continuum result. 
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Figure 2: A comparison of the quadratic potential for different blocking prescriptions in a cubic 
box of size 16 3 . A non-physical constant was subtracted to match the curves. 



2 We have checked this explicitly by constructing FP Polyakov loops in the linear approximation. 
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3.2 The temperature scale 



In this section wc determine the critical couplings (3 C for the colour deconfinement phase transition 
on lattices with different temporal extensions N t . As in ref. Q] we use the definition of the Columbia 
group |50|: For a given lattice size N% x N t we perform MC simulations for various couplings /3 
around the expected transition region measuring the Polyakov loop averaged over the lattice: 

-, Nt-l 

P =^E Tr II Mn^ft) (3.4) 

s n n =0 

We evaluate the angle 9 from tan(#) = Im(P)/Re(P) and measure the fraction of the time /20 
that 9 lies in the range of ±20° of the directions of the Z(3) roots. Finally, the phase transition is 
defined as the coupling /3 where the simulations yield for the deconfinement fraction 

3 1 

fd = g/ao - 2 ( 3 - 5 ) 

the value 0.5. 



We measured the critical couplings for lattices of temporal extensions Nt — 2, 3, 4 and 6 
with different spatial sizes N s . The results are listed in table [| The deconfinement fraction has 
been determined at various couplings in the vicinity of the phase transition. At each coupling 
we have performed from 2000 to 30000 sweeps. The critical coupling j3 c was obtained by a linear 
fit of the fractions bracketing fd{Pc) = 0.5. Its error was obtained by a jack-knife analysis. The 
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3.361(5) 
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3.378(3) 
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3.385(9) 
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3.395(3) 




3.678(3) 
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3.399(5) 


3.581(4) 
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3.686(3) 


3.91(4) 


12 3 




3.587(5) 


3.691(5) 


3.87(1) 


143 






3.691(5) 


3.882(7) 


16 3 








3.882(8) 


00 


3.400(3) 


3.588(4) 


3.695(4) 


3.886(13) 



Table 5: Critical couplings at finite volume and extrapolated to infinite volume for the FP action 
with parameters Ilia in Table ||. 

critical couplings in the thermodynamic limit (A^ s — *■ 00) were obtained by extrapolating the finite 
couplings using the finite size law fl£ , 30 1: 

3 



p c (N t ,N s ) = (3 c {N t , 00) - c(JVt) 
The finite size scaling behaviour is plotted in fig. [}| 



(3.6) 



3.3 The full potential 



To repeat similar measurements of ref. we have measured the static gg-potential V(r; T) at finite 
temperature T = |T C for the type III action. The result for j3 — f3 c (N t = 2) = 3.40 on a lattice 
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Figure 3: Finite size scaling for the critical couplings of the deconfinement phase transition using 
the parametrization Ilia. The j3 = oo extrapolations are included. 



6 3 x 3 is shown on fig. ||. The constant in the potential is fixed by setting V(r = 1/T C ;T) = 0. 
For the Wilson action we used 12 3 x 6 and 6 3 x 3 lattices, again at T = |T C . These data and 
those for type I action are taken from rcf. [Q. As opposed to the Wilson action, both FP actions 
for Nt = 3 have a proper slope at large distances. One also observes that the type III action has 
a somewhat smoother potential than type I. The difference is, however, not as striking as for the 
perturbative potential on fig. |^. Note also a phenomenon which looks strange at first sight. The 
third point from the right on the type I and type III plots — which corresponds to the diagonal 
distance (2,2,2) — has larger error bars than the others. In the Wilson case this point is not plotted 
because the measured correlation value was lost in the noise. This is also a sign of strong violation 
of rotational symmetry, and is not completely restored in the parametrized FP actions either. The 
reason for this phenomenon (besides the trivial factor of different multiplicities) is, perhaps, that 
in the Wilson action there is no direct interaction between diagonally separated links, and in the 
actual parametrization of the FP actions used here no spatial diagonal interaction term is included 
- this would be an 8-link 'parallelogram'. (Such direct interaction term is certainly present in 
the true FP action.) As a result, the diagonally separated Polyakov loops are not bound to each 
other so strongly as they should be. A similar behaviour is observed for the distance (1,1,1), but 
the errors are too small there to be seen in the figure. 
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Figure 4: The gg-potential V(r;T)/T c vs. rT c for different actions at T = 2/3 T c . Upper left: 
Wilson, upper right: type I, lower left: type III action. These three plots correspond to N t = 3. 
The dashed line shows the function — 0.24/r + 2.48r — 2.19 , and serves to guide the eye. The lower 
right graph shows again the potential for type III and Wilson action for N t — 3, together with the 
N t — 6 result for the Wilson action. The spatial size in all cases is N s — 2N t . 
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4 Conclusion 



Using a more general smearing kernel we obtained a new FP action which in the quadratic approx- 
imation performs better than those obtained previously — it has a somewhat faster decay rate 
for the interaction coefficients and produces a perturbative potential close to the continuum one. 
As far as the cut-off effects for finite (3 are concerned, the performance of this new FP action, or 
rather of the actual parametrization suggested here, has to be investigated further. To improve 
the parametrization it could also be useful to include rough classical solutions (like instantons or 
constant Abelian fields) with analytically known action values fL^,|l4||. 

As mentioned in the introduction, A. Papa has measured the cut-off effects in the free energy 
for type I and type Ufa actions ||(| . For the type Ilia action at T = 2T C on a lattice with temporal 
extension N t — 3 no cut-off effect has been observed within the small errors — the result agrees 
with the continuum prediction of ref. Jl^,|l9| and even at N t = 2 the error is ~ 10%. For the type I 
the deviation is ~ 10% at Nt — 3. This has to be compared to the Wilson case where the cut-off 
effect is ~ 100% at Nt = 2 and even at Nt = 4 it is still ~ 20%. Because the necessary simulation 
time grows very fast with N t — about -/V t 10 — it is essential to keep the cut-off effects as small as 
possible. 
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